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Abstract 

In supersymmetric scenarios with a long-lived stau, the LHC experiments 
provide us with a great environment for precise mass measurements of superpar- 
ticles. We study a case in which the mass differences between the lightest stau 
and other sleptons are ~ 10 GeV or larger, so that the decay products of heavier 
sleptons are hard enough to be detected. We demonstrate that the masses of 
neutralinos, sleptons, and squarks can be measured with a good accuracy. 



1 Introduction 



Signatures of supersymmetry (SUSY) at the LHC experiments crucially depend on 
what the lightest supersymmetric particle (LSP) is. Many of studies have assumed 
a neutralino as the LSP motivated by a possible explanation of dark matter of the 
universe. In this case, final states of SUSY events will be accompanied by missing 
momentum carried away by two neutralinos. The mass determinations of SUSY par- 
ticles in such cases are not a straightforward task. One needs to combine various 
measurements to extract the masses out of observables. 

Situation drastically changes when we assume that the scalar tau lepton (stau) 
is lighter than the lightest neutralino and thus long-lived. First, one can discover 
the stau by looking for anomalous tracks at the inner tracker and the muon detec- 
tor [TJ 121 [3] . Their momenta and velocities can be measured by analyzing the tracks, 
with which one can measure the stau mass with a good accuracy jU El El [7j. The 
momentum information enables us to perform precise determination of the properties 
of SUSY particles such as masses of other superparticles [HJ (H [TDJ, [TTJ , the spin of the 
stau H2], the lifetime of the stau [Hfl QU EH [TB] and P/CP/T violations in the SUSY 
interactions [T7] . 

The long-lived stau is not just motivated as a golden scenario for collider exper- 
iments. There are many underlying models and parameter spaces in them to realize 
the scenario. Examples include the familiar ones such as supergravity and gauge me- 
diation models. The lifetime of stau depends on details of the microscopic models; 
the stau can decay into a gravitino and a tau lepton if kinematically allowed or into 
standard model particles if i?-parity is violated. There is a cosmological constraint 
on the stau lifetime from Big-Bang nucleosynthesis, but it can easily be evaded unless 
the lifetime is extremely long such as greater than 0(1000 sec) [18J. 

Once the stau tracks are discovered and the measurement of the stau mass is done 
at the LHC, the next important quantities to measure will be other superparticle 
masses. In particular, measurements of selectron and smuon masses tell us about 
how special the third generation is in the microscopic interactions between the SUSY 
breaking and the standard model sectors. Even if the strength of the Yukawa inter- 
action is the only difference between the stau and the first two generations, the mass 
splitting contains the information on the size of the Yukawa interaction (i.e., tan/3) 
and also the size of the quantum corrections to the masses (i.e., the messenger scale). 
An even more interesting quantity will be the mass splitting between the selectron 
and the smuon. Since the Yukawa coupling constants are small, a large splitting of 
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more than 0(100 MeV) level requires an explanation from a microscopic theory. 

The mass measurements of superparticles in the long-lived stau scenario was first 
studied in [8]. They used a parameter point in the gauge-mediation model where the 
selectron, the smuon, and the stau are all long-lived due to a small mass splittings 
among three sleptons (nig ~ ~ m f ). Recently, the case with a long-lived selec- 
tron with a nearly degenerate smuon (m g ~ < nif) was studied in there it 
was proposed to measure the mass splitting between the smuon and the selectron by 
looking at the peak locations of the ee and fie invariant masses. 

In this paper, we propose a method to measure the masses of superparticles in the 
case where m f < m e - ~ < m x o in the long-lived stau scenario at the LHC. When 
the mass differences between the stau and other sleptons are larger than ~ 10 GeV, 
the decay products are hard enough to be detected, contrasted to the previous studies. 
This level of mass splitting is a natural expectation in many models especially for large 
values of the tan f3 parameter (which is somewhat motivated by the Higgs boson mass 
bound). 

At the LHC, staus are mainly produced at the last step of the cascade-decay chains 
such as the neutralino decay x° ~~ > T ^ an d the chargino decay X s1 tv. It has been 
demonstrated in various models that the neutralino mass(es) can be measured by 
looking for an endpoint of the j T -f invariant mass distribution [SJ El EI] (where j T is 
r-tagged jet). We first follow this analysis and point out usefulness of the charge in- 
formation of the r-jets in this study. We next show that the measured neutralino mass 
in turn can be used to determine the smuon and selectron masses by reconstructing 
the decay chains: x° — > ^ ~~ * A t (^ r / i ) an d X° ~ > e & ~ ¥ e(fre) followed by a r decay. 
We use the hadronic r decay for the analysis. Although there are invisible neutrinos 
from r decays, four-momentum of r can be reconstructed in the event-by-event basis 
by using the knowledge of the neutralino mass under an assumption that the neutrino 
is emitted along the direction of the r-jet (which is valid when r is highly boosted). 
The smuon and selectron masses can then be measured in each event up to a combi- 
natorics. Therefore, we can see sharp peaks at the smuon and selectron masses in the 
distributions of the frfi and f re invariant masses, respectively. We demonstrate that 
the masses and the mass difference can be measured with accuracies of 0(100 MeV). 
This level of accurate mass measurements will provide us with a very important infor- 
mation on the underlying microscopic theory. We also proceed to reconstruct squark 
masses by using the q — > q\\ — > qff decay. We can see a sharp peak in the q-r-f 
invariant mass distribution. 

The organization of this paper is as follows. In Section [21 basic ideas to measure 
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the neutralino, slepton and squark masses are explained. In order to demonstrate that 
our basic ideas work, we show the results of Monte Carlo (MC) analysis in Section [3j 
We will discuss implications to microscopic theories in Section HI Section [5] is devoted 
to conclusions. 

2 Basic Ideas 

As mentioned in Introduction, in a class of SUSY models, the lightest stau fi can be 
the lightest superparticle in the MSSM sector (which we call the MSSM-LSP). The 
lifetime of f\ can be long enough to escape the detector before their decays, thereby 
one can treat it as a stable particle in collider experiments. 

If fi is long-lived, we expect a unique LHC phenomenology, very different from 
the case of the neutralino LSP. In particular, we will see tracks of f\ and thus the 
momentum information on the MSSM-LSP will be available in contrast to the missing 
momentum in the neutralino LSP case. The momentum information enables us to 
easily reconstruct SUSY events. It is also notable that, once the fi track is identified 
in an event by looking at both the velocity and the momentum, the event can be 
distinguished from the standard-model events, which in principle eliminates standard- 
model backgrounds. It has been studied that the selection of slow tracks effectively 
reduce the muon background [6j [7]. In the following analysis, we assume that the f\ 
track can be distinguished from the muon track if (3^ < 0.9 [6], and we neglect the 
standard-model background. 

Collider phenomenology of the long-lived f\ scenario will be quite different depend- 
ing upon the mass spectrum of the superparticles. In this paper, we consider the case 
with the following mass relation: 

where m? 1 , rng R , m^, and m x o are masses of fx, lighter selectron Br, lighter smuon 
[Mr, and the lightest neutralino Xi-> respectively. In addition, we assume the following 
for simplicity: 

• All the colored SUSY particles are heavier than Xv 

• The lighter sleptons are (almost) right-handed. 

• The heavier sleptons, which are almost left-handed, are heavier than Xv 

We pay particular attention to the case that the mass differences rag R ^ R — m fl and 
m x o —rrie Rt fi R are both sizable; if so, the decay products of Xi an d sleptons are energetic 
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enough to be detected. These conditions are realized in a wide class of SUSY breaking 
models; one of the examples is the minimal gauge-mediated SUSY breaking (GMSB) 
model [23] with a large value of tan f3. 

With the mass spectrum mentioned above, some of the neutralinos (in particular, 
Xx) decay as x°i ~^ and x°i ~^ (Here and hereafter, I stands for e and /i.) 

In the latter case, Ir then decays into three-body final state: Ir — > It ff. Since the 
momentum of f\ can be measured, we can reconstruct mg R , irijx R , and m x o using these 
decay processes. 

The basic procedures are as follows. For the reconstruction of the neutralino mass 
m x o, one can use the decay process Xi ~ * T T i, followed by the hadronic decay of r. 
Since r-jets have a distinguishable feature from the typical QCD jets; a very narrow jet 
containing small number of charged track(s), one can identify it with a high efficiency 
(Here, we only use 1- and 3-prong decay of r.) If we consider the distribution of an 
invariant mass: Mj T ? = \/ (pj T + Pf ) 2 , where pj T and Pf are four-momenta of r-jet and 
f, respectively, there should be an upper endpoint at m v o. 

Once the lightest neutralino mass is measured, it can be used to reconstruct the 
slepton masses mg R and For this purpose, we use the decay process Xi ~~ > 

(with I = e or /i) followed by — > / =F rf 1 and r —> j T v T . Since the r leptons from the 
In decay are expected to be highly boosted, the directions of r and r-jet are (almost) 
aligned. Then, the four momentum of r is obtained by requiring (pi++pi- +p T +Pf) 2 = 
m 2 . Once p T is known, we define Mj± = \J {pf + p T + p f ) 2 . Because one of the 
charged leptons is from the decay of Z^, M~ l+ or is equal to rni R - Therefore, we 
expect to have a sharp peak at the slepton masses in the distribution of M~ l± . 

In order to reconstruct the neutralino and slepton masses with the above-mentioned 
procedure, we need enough amount of neutralinos to perform the statistical analysis. 
The production cross section of neutralinos, of course, depends on the model param- 
eters, but in many cases the neutralinos are copiously produced from the decay of 
squarks. In fact, squark mass measurements may also be possible by reconstructing 
the decay products as we demonstrate later. 

In Ref. [8], the squark mass measurement was discussed for the case that Br and 
£lr, as well as f 1; are long-lived. In such a case, the momentum of Ir can be well 
measured. By using the decay chain q — > qxi — > qllR, it was pointed out that a 
selection based on the invariant mass of the system {I^Jr 1 ) works to find / and Ir 
which are from the same neutralino. Then, the invariant mass of the system (j, I, Ir) 
(with j being one of high-p^ jet) directly provides the squark mass. 

If Ir decays in the detector, this method is not applicable. Even in such a case, 
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however, the squark mass information can be obtained by reconstructing the momen- 
tum of r. Considering the lightest neutralino production process q — > q\\ — > qrf\ 
followed by hadronic decay of r, and using the fact that m x o can be determined from 
the endpoint analysis, we can first reconstruct the momentum of r by requiring that 
the invariant mass a/ (p T + Vf ) 2 be equal to rn x o. Then, the squark mass can be ob- 
tained by a/ (Pq + Vt + Pf) 2 , where p q is the momentum of the quark from the squark 
decay, whose information is in principle imprinted in the momentum of one of the ob- 
served jets. We note here that the flavor information on the primary quark is hardly 
obtained except for the 6-jet, so we can perform only the flavor-blind analysis. In a 
large class of models, however, first- and second-generation squarks are almost degen- 
erate. (For example, in the GMSB and mSUGRA models, this is the case.) In such a 
case, as we will see in the next section, a clear peak corresponding to the squark mass 
is obtained even though we cannot specify the flavor of the jets. 

Of course, in the actual situation, the measurements of the masses are not straight- 
forward. This is because (i) the momenta of the decay products of Xi are measured 
with some uncertainties, (ii) jets with small multiplicity (i.e., r-jet-like objects) are 
also produced by the QCD process, which mimics the hadronically decaying r-lepton, 
and (iii) combinatorial backgrounds should exist. In the next section, we discuss how 
well our idea works at the LHC experiment using MC analysis, taking account of the 
effects of (i) — (iii). 

3 Numerical Analysis 

We demonstrate in the following the method to measure the sparticle masses presented 
in the previous section by performing a Monte Carlo simulation. As a model with the 
spectrum in Eq. (I2.ip . we use the minimal GMSB model. The superparticle spectrum 
is parametrized by A (the ratio of the F- and A-components of the SUSY breaking 
field), M mcss (the messenger scale), N 5 (the number of messenger multiplet in units of 
5 + 5 representation), tan/? (the ratio of the vacuum expectation values of two Higgs 
bosons), and the sign of /i (the SUSY invariant Higgs mass). We take 

A = 60TeV, M mess = 900 TeV, N 5 = 3, tan/3 = 35, sign(/i) = +. (3.1) 

Even though we adopt specific models for our analysis, it should be noted that our 
procedure works in a wider class of models as far as the mass relation (12.11) holds. 

The mass spectrum of superparticles is calculated by using ISAJET 7.64 [19]. In 
order to see how effective our method is, we reduce the mass of fx,R by 1 GeV. With the 
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present choices of parameters, the MSSM-LSP is f±, and the lightest neutralino, which 
is almost the Bino, is heavier than lighter sleptons cr and (Ir. The mass spectrum of 
the superparticle is summarized in Table [D The branching ratios of the decay of Xx 
are given by Br(xx ~~ > ttx) = 61.4 % and Br(xx — > 11 r) = 19.3 %, while Ir dominantly 
decays as l R — > Itt\ (~ 100 %). 

We have generated 66960 SUSY events corresponding to the luminosity of 100 fb _1 
at a pp collider with the center of mass energy of 14 TeV by using the HERWIG 6.510 
package [20l [21]. (The total cross section of the SUSY events is 669.6 fb.) Events are 



passed through the PGS4 detector simulator [22] cj The momentum resolution of fx 



is assumed to be the same as those of muons as long as 0.4 < (3^ < 0.9. 
3.1 Neutralino masses 

We first discuss the neutralino mass measurement explained in the previous section. 
We identify a decay process: 

X? - 7*7?, (3.2) 



followed by the hadronic decay of the r-leptonr 2 I The following selection cuts are 
applied: 

la) At least one f% with the velocity 0.4 < (3 fl < 0.9; such f% is assumed to be 
detected with the efficiency of 100 % with no standard-model background. In 
the study of the invariant-mass distribution, fx with the velocity in this range 
are used. 

lb) At least one r-tagged jet with p T > 15 GeV, which is denoted as j T . 

By using events passed the selection cuts, we calculate the invariant mass for all the 
possible pairs of (j T ,fx)' 



M jTfl = ^(p jT +Pr) 2 . (3.3) 

The charges of j T and fx, which are both observable, should be opposite for signal 
events; we call such events as opposite-sign (OS) event. In Fig. [TJ (left), we plot the 



^In PGS, effects of energy leakage into the hadronic calorimeter are included for electromagnetic 
objects (i.e., and 7), which results in an underestimation of the energy of e ± . We expect that 
the energy of will be calibrated by requiring that the Z-boson mass is well reconstructed. So, we 
estimate the observed energy of by neglecting the energy leakage. 

# 2 We have also considered the possibility to use the leptonic decay mode of r-lepton. However, for 
the signals from the leptonic decay mode, the combinatorial background is so severe that the results 
are much worse than the case with the hadronic decay mode of r. 
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Particle Mass (GeV) 
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1234.28 
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1180.19 


ti 


1082.85 


h 


1195.08 




1145.24 


h 


1185.83 




388.05 


h 


396.19 


h 


402.57 




383.80 


en 


194.39 




193.39 


T\ 


148.83 


x\ 


239.52 


xl 


425.92 


xl 


508.41 


xl 


548.67 


xf 


425.45 


xt 


548.43 


h 


115.01 



Table 1: Masses of the superparticles and the lightest Higgs boson h in units of GeV. 
The input parameters are A = 60 TeV, M mess = 900 TeV, N 5 = 3, tan/3 = 35, 
sign(/i) = +. (We use the top-quark mass of 171.3 GeV.) We reduced the mass of [Lr 
by 1 GeV. 
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Invariant Mass (stau tau) 



Fitting : M 



bino 
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X 2 / ndf 28.779 / 46 
A 735.84 ± 257.06 
M 241.11 ±0.96203 
a 3.4819 ± 1.32 
N 22.662 ± 5.3327 




260 
mv [GeV] 



Figure 1: [Left]: The invariant mass distributions of opposite sign (black) and same sign 
(red) pairs for the GMSB case. Here, we take C = 100 fb _1 , and the postulated stau mass 
is taken to be the underlying value. [Right]: The fitting of the lightest neutralino mass after 
the charge subtraction. 



distribution of the invariant mass Mj rfl . One can find a sharp drop-off at M fl j T ~ 
240 GeV, which is close to the input value of the lightest neutralino mass. However, 
one can also see a long tail above the drop-off. Those backgrounds come from fake 
r-jets (mis- identified QCD jets) as well as from wrong combination where r and f 
have different parents. Since those backgrounds are charge-blind, their contributions 
to the OS and same-sign (SS) histograms are expected to be the same amount. In 
Fig. [I] (left), we also show the distribution of the invariant mass for the SS event. One 
can see that the number of OS events is significantly larger than that of SS events for 
M fl j T < 240 GeV, and those two become comparable for a larger invariant mass. This 
fact confirms our expectation that the tail is mostly from the fake r-jets and wrong 
combination. 

By taking a difference between the OS and SS events, one can subtract the contri- 
butions from the backgrounds. The distribution of the signal events after the charge 
subtraction is shown in Fig. [T] (right). We can see a clearer endpoint. We search 
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Invariant Mass (stau tau) 



Fitting : M 
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M inv [GeV] 



Z 2 /ndf 20.415/36 
A 31 7.75 ±146.54 
M 430 ± 4.4684 
a 9.8908 ± 4.8052 
N 15.504 ±3.8853 




350 400 



450 500 
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Figure 2: [Left]: The invariant mass distributions of opposite sign (black) and same sign 
(red) pairs for the GMSB case. Here, we take C = 100 fb _1 , and the width of the bin is 
10 GeV. [Right]: The fitting of the second-lightest neutralino mass after charge subtraction. 



the upper endpoint by fitting the histogram with the Gaussian-smeared triangular 
function (including the effect of background) [24J : 



(Number of events) = A I dz exp 



+ N BG , (3.4) 



where A, M^ max \ a, and Nbg are fitting parameters; in particular, M*™ 3 ^ corresponds 
to the upper endpoint. The endpoint is determined to be (241.11 ± 0.96) GeV, which 
is consistent with the underlying value of m x o. 

In the above study, the charge subtraction method was very powerful in finding 
the neutralino endpoint. One can further try to find the second-lightest neutralino, 
which is almost the neutral Wino in this model. Although no clear endpoint can be 
seen in the OS events at ~ m x %i if shows up after the subtraction. (See Fig. [2], where 
the distribution of Mj Tfl is shown with the width of the bin of 5 GeV.) By fitting 
the edge with the same function, m x o is measured to be (430.0 ± 4.5) GeV, which is 
slightly larger than the underlying value. 
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Before closing this subsection, we comment on the -dependence of the recon- 
structed mass. The stau mass is expected to be measured in the long-lived stau 
scenario by combining the velocity and the momentum informations on the f\ track; 
the expected error in the stau mass measurement is Sm^ ~ 0(100 MeV) [61 [7J. The 
neutralino mass measurements will therefore be affected by 0(100 MeV), which is less 
significant compared to the statistical uncertainties. 



3.2 Slepton masses 

We proceed to the discussion of the slepton mass measurements. In particular, it is 
interesting to measure the masses of e# and fiR since one can learn the flavor structure 
of the model. Note that study of those particles are difficult in the Bino LSP scenario 
since they do not appear in the cascade decays. In order to measure mj , we use the 
decay chain 

x l -» -> PFrfu (3.5) 

followed by the hadronic decay of r. In the decay chain (I3.5p . the charges of r and f 
are opposite. We apply the following selections: 

2a) At least one f\ with the velocity 0.4 < < 0.9. 

2b) At least one r-tagged jet with pr > 15 GeV. 

2c) At least one pair of isolated opposite-charge same-flavor leptons. We require 
Pt > 15 GeV for leptons. 

For each event, we choose all the possible combinations of l~,j T , ^i) with requiring 
charges of f\ and j T to be opposite. Assuming that l + , l~, j T , and f\ are all from 
the lightest neutralino whose mass is determined already, we can calculate the four- 
momentum of the r-lepton. Because the r-lepton from the slepton decay is usually 
ultra-relativistic, r and j T are expected to be emitted to (almost) the same direction. 
The four-momentum of r can thus be estimated as 

p T = z~ 1 Pj T , (3-6) 

where 

2 Pj T ■ (Pi+ + Pi- + Pf) 



mlo ~ (Pi+ + Pi- + Pf) 2 ' 



(3.7) 



with fh x o being the postulated value of the lightest neutralino mass for the analysis. 
The events with z > 1 are rejected. In order to reconstruct the slepton mass wif , we 
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Figure 3: Distribution of Mi with I = e (left) and \x (right) for the GMSB case with 
£ = 100 fb _ . Here, the postulated stau and neutralino masses are chosen to be the 
underlying values. 



study the distribution of the following quantity: 

Mj = ^(pi+Pr + z-^) 2 - (3.8) 

Here, pi is the four momentum of one of two charged leptons in (/ + , l~,j T , fx). Because 
we cannot tell which lepton is the one from the slepton decay, we calculate the invariant 
masses by using both of the possibilities. 

In Fig. [31 we show the distributions of with I = e and /1. Here, the postulated 
values of the stau and the lightest neutralino masses are taken to be equal to the un- 
derlying values. Even though the distribution contains the combinatorial background, 
one can see sharp peaks at the slepton masses. Fitting the peaks with the Gaussian 
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function plus a constant background, 



(Number of events) = A exp 



L(Mj-M^) 2 +N BG , 



(3.9) 



with A, M^ eak \ a, and Nbg being parameters, the peak positions j\/f( peak ) are deter- 
mined to be (194.46±0.11) GeV and (193.44±0.07) GeV for I = e and ft, respectively. 
Those values are in very good agreements with the actual slepton masses. 

In the slepton mass determination, we should consider errors associated with the 
uncertainties in the mass measurements of stau and the lightest neutralino. In partic- 
ular, the reconstructed slepton masses are sensitive to the uncertainty in the lightest 
neutralino mass. In Fig. HI we show the reconstructed slepton masses (i.e., the l-a 
upper and lower bounds with C = 100 ftT 1 ) as functions of the postulated lightest 
neutralino mass. One can see that the reconstructed slepton masses depend linearly 
on the postulated lightest neutralino mass. As mentioned in the previous section, the 
lightest neutralino mass can be measured by the endpoint study, and the uncertainty 
is ~ 1 GeV, which gives the error of [Smj }sm ~ 1 GeV in the slepton mass mea- 
surement. Importantly, however, the mass difference mg R — ttj^ is insensitive to the 
neutralino mass. The dependence of the reconstructed slepton masses on the pos- 
tulated stau mass is proportional to (ui^o — m^)/m x o, and thus rather weak in the 
present model. 

3.3 Squark mass 

Sqark masses can be also measured by reconstructing the decay chain: 



followed by hadronic decay of r. In order to use the decay chain (13.101) . we adopt the 
following event selections: 

3a) At least one fx with the velocity 0.4 < < 0.9. 
3b) At least one jet with px > 100 GeV. 
3c) At least one r-tagged jet with px > 15 GeV. 
3d) No isolated lepton with px > 15 GeV. 

The condition 3d) is to eliminate the mis-reconstruction of leptonically decaying r. 

Then, we consider all the possible combinations (j,j T ,Ti), where j is one of four 
leading jets with px > 100 GeV. Assuming that j T and f% are both from the decay of 




qrr, 



(3.10) 
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Figure 4: The reconstructed slepton masses as functions of the postulated lightest 
neutralino mass for the GMSB case. 

the neutralino, we first reconstruct the momentum of r, which is given by 

Pr = Zq 1 Pj T , (3.H) 

where 

H = 2 2 JT ' P \ - (3-12) 

TOo — TO- 

Combinations with Zq > 1 is eliminated. Then, we study the distribution of the 
following variable: 



M^y/to+Pf + zfp^)*. (3.13) 

The distribution of is shown in Figs. [5] for OS and SS events. We can see 
a clear peak for the OS event at the position corresponding to the squark masses, 
while no clear peak is seen for the SS event. By fitting the histogram after the charge 
subtraction with the Gaussian function, the peak position is obtained to be trig = 
(1172.1 ± 1.2) GeV, while underlying squark masses are mQ R = 1184 GeV, mj 
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Figure 5: [Left]: The invariant mass distributions of opposite sign (black) and same sign 
(red) events. Here, we take C = 100 fb _1 , and the postulated stau mass is taken to be 
the underlying value of the stau mass. [Right]: The fit of the number of events with the 
Gaussian function. 



1180 GeV, m 



"L 



1232 GeV, and mg L 



1234 GeV. The dominant contribution to 



the histogram is from right-handed squarks. Compared to the right-handed squark 
masses, the observed peak location is smaller than the underlying values by ~ 10 GeV. 
It is partly from the mis-measurement of the jet energy and also from effects of the 
background. The uncertainties from the mis-measurements of the Bino and the stau 
masses are much smaller than the statistical error as far as the luminosity of C ~ 
100 fb _1 is used. 

In principle, a similar analysis can be performed by using the second-lightest neu- 
tralino (which is almost the Wino), which should provide information on the left- 
handed squark mass. However, the left-handed squarks decay into chargino (plus 
jet) as well as into neutralino. In addition, it is difficult to distinguish the left- and 
right-handed squark production events in the event-by-event basis. These become the 
source of extra background and reduce the signal-to-background ratio. Thus, for the 
study of the left-handed squark mass, we should perform more sophisticated analysis. 
A detailed discussion will be given elsewhere [25] . 
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4 Implications 



We have seen that, in the long-lived stau scenario, masses of various SUSY particles 
are expected to be measured with a very good accuracy at the LHC. This level of 
precise deterimnation of the mass spectrum may reveal the origin of SUSY breaking 
terms. 

As we have seen, the masses of the lightest and the second-lightest neutralino, 
which almost correspond to the Bino and the neutral Wino, respectively, can be deter- 
mined with the endpoint analyses. Then, combined with the gluino-mass information 
from, for example, the cross-section information for the process pp — > gg, we can test 
the GUT relation among the gaugino masses. 

The masses of cr and JIr are also measured with good accuracies. Even though 
the mass determinations of cr and flu are sensitive to the uncertainty in the lightest 
neutralino mass, the measured value of the mass differene mg fl — m^ fl is not affected 
to it. The mass difference mg fl — mp iR contains various informations. Neglecting the 
flavor mixing, the slepton mass matrix has the form: 

M\4 ml l^ A «« ), (4.14) 

V ^l,LR m R + ^l,RR J 

where mi and rriR are independent of the flavor index. Then, assuming that A <C 
m 2 LRl the mass difference is given by 

m " - m >« H i (A « R - - 2m a (4-^) (Ag 4" - A *4«>- < 4 ' 15 > 

Thus, the slepton-mass difference is sensitive to the left-right mixing, whose infor- 
mation is in Aj LR , as well as to the difference of the diagonal elements of the mass 
matrix. Even if the dominant contribution to the masses of sleptons (with the same 
gauge quantum numbers) is almost universal at some high energy scale, the univer- 
sality is affected by various effects. 

For example, in the GMSB models, the non-universality is originated from the 
supergravity effects and also from one-loop corrections with the muon Yukawa inter- 
action. The supergravity effects are estimated as Aj LL ~ Aj RR ~ 0(m^ 2 ), with m 3 / 2 
being the gravitino mass. Thus, if the gravitino mass is as large as a few GeV, which is 
merginally consistent with the flavor-violation constraints depending on other SUSY 
parameters, the supergravity effect may be seen assuming that the mass difference 
rrie R — mp iR is determined with the accuracy of 0(100 MeV). 

The size of the loop effects is enhanced when tan f3 is large, and the contribution 
to the mass splitting is roughly estimated to be 0(100 MeV) for tan/? > 30. The 
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Figure 6: The selectron-smuon mass splitting calculated in the sweet-spot model. 

actual value of the loop effects depends on the mechanism to generate the Higgs mass 
parameters since there is a one-loop diagram with Higgs fields in the loop. As an 
example of explicit models we calculate the mass splitting by solving renormalization 
group (RG) equations in the sweet-spot supersymmetry model [10]. SUSY breaking 
terms are determined by the parameters: A, M mess , A 5 , and \i. The Higgs mass 
parameters are generated at the unification scale rather than the messenger scale in 
this model, and hence there is a large logarithm in the quantum corrections. Also, 
the tan/3 parameter is predicted to be large due to the boundary condition at the 
messenger scale, B = 0. The Yukawa coupling constants for the charged leptons are 
therfore enhanced. For these reasons, the mass splitting, Am = — m^ fl , is expected 
to be large. Note also that for the same reason, f\ can be significantly lighter than 
the other sleptons, that enables us to perform the analysis in the previous section. 
The mass spectrum is calculated with fixing A and N 5 as 

A = 60 TeV, A 5 = 3, (4.16) 

so that the gaugino masses are similar to the ones in the example we took in the 
previous section. Since the value of /i (the Higgsino mass parameter) is strongly 
correlated to the stau mass in this model (through the RG effects mentioned above), 
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we can trade the input parameter /i with . We show in Fig. [6] the mass splitting 
as a function of the stau mass for various values of the messenger scale. For M mess > 
10 10 GeV, the slepton is heavier than the Bino with this set of parameters. The lines 
are terminated at the point where the stau becomes heavier than the other sleptons. 
One can see that the mass splitting can be as large as ~ 300 MeV which is large enough 
to be observed at the LHC This calculation demonstrates that the mass splitting has 
a rich information on the microscopic theory to generate the SUSY breaking terms, 
especially those for the Higgs fields. It is also interesting to note that the LHC may 
be able to see an effect of the muon Yukawa coupling. 



5 Conclusions 

We have developed methods to measure the superparticle masses at the LHC in the 
long-lived stau scenario. We have concentrated on the scenario where m fl < mj R < 
m x o (with I = e and fi) and demonstrated that the masses of neutralinos, sleptons 
and squarks can be well determined by endpoint or peak analysis. 

In the neutralino mass measurements by the endpoint analysis of the Mj T f invariant 
masses, we have seen that the charge subtraction method is useful to identify the 
endpoints. In the sample point we have adopted, the estimated error in the lightest 
neutralino mass measurement is ~ 1 GeV. We have shown that even the second- 
lightest neutralino mass can be measured with an accuracy of 4 — 5 GeV. 

Once the lightest neutralino mass is known, it can be used to determine the masses 
of selectron and smuon. With the decay chain Xi —> I ^ ~~ > / ± / =F rf 1 , followed by 
hadronic decay of r, we first reconstruct the momentum of r by using the relation 
m x o = a/ (pi+ + pi- + p T + Pf ) 2 - Then, the slepton mass mf is determined in each 
event up to a combinatorics. With this method, we have seen that very sharp peaks 
are obtained around the underlying values of the slepton masses, which gives precise 
determination of the slepton masses. We have estimated the error in the slepton mass 
determination, which is ~ 100 MeV. Such a precise measurement of the slepton mass 
enables us to study effects of renormalization group, supergravity, and/or left-right 
mixing on the slepton masses. We also demonstrated that a sharp peak corresponding 
to (right-handed) squarks can be observed by using the q — > qxi — ► qTf\ events. 
Acknowledgments: We thank Prof. S. Asai for useful discussion. This work was sup- 
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cation, Science, Sports, and Culture of Japan, no. 21840006 (R.K.) and no. 19540255 
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